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Abstract 



'^ ' The influence of long-range Coulomb interactions on the properties of one- 

^ ' dimensional (ID) strongly correlated electron systems in vicinity of the metal- 

insulator phase transition is considered. It is shown that unscreened repulsive 

}^ • Coulomb forces lead to the formation of a ID Wigner crystal in the metallic 

C^ \ 

phase and to the transformation of the square-root singularity of the compress- 
ibility (characterizing the commensurate-incommensurate transition) to a log- 
arithmic singularity. The properties of the insulating (Mott) phase depend on 
the character of the short- wavelength screening of the Coulomb forces. For a 
sufficiently short screening length the characteristics of the charge excitations 
in the insulating phase are totally determined by the Coulomb interaction 
and these quasipartic les can be described as quasiclassical Coulomb solitons. 



The metal-insulator transition induced by strong correlations in an electron system is a 
problem of permanent theoretical interest. In recent years, this problem became particu- 
larly important in connection with the discovery of high-temperature superconductivity and 
subsequent attempts to develop a consistent theory of this phenomenon [1] . 

Much progress has been achieved in one-dimensional models where exact results or well 
developed non-perturbative methods are available, making it possible to obtain an analytical 
description of the dynamics of the Mott-Hubbard transition in a wide range of bare parame- 
ters of the system. For example, using the Bethe-ansatz solution of the ID Hubbard model, 
an analytical expression has been recently obtained for the charge stiffness of a finite-size 
system close to the transition as a function of the on-site repulsion U [2] . 

It is well known that, at half filling (one electron per site), the ID repulsive Hubbard 
model describes a Mott insulator. In this case, the charge excitations have a gap in their 
spectrum and at weak Hubbard interaction, U <^t {t being the bandwidth), can be regarded 
as topological solitons (kinks) of the quantum sine-Gordon model with coupling constant 
13"^ = 8tt (see e.g. [3]). At finite deviations from half filling the system is in a metallic phase, 
but the charge carriers ("holons") keep the memory of the Mott phase; they can still be 
considered as solitons with a characteristic size proportional to the correlation length C,o in 
the insulating phase. At low densities, nC,o <^ 1, the holons interact weakly and can be 
described as free massive spinless fermions [2,3]. 

In this paper we study the effects of a long-range Coulomb interaction on the metal- 
insulator transition in strongly correlated one-dimensional electron systems. We shall as- 
sume that the Mott phase for the charge excitations is described in terms of the sine-Gordon 
model, while the transition to the metallic phase is related to creation of solitons at the chem- 
ical potential exceeding the charge gap. Near the phase transition, interactions between the 
solitons can be neglected, so that in the absence of Coulomb correlations the properties 
of the system near the transition are completely determined by Fermi statistics of quan- 
tum solitons. In particular, the compressibility of the system has a trivial (in fermionic 
language) square-root singularity typical for the commensurate-incommensurate transition 



[4,5,6]. Notice that, within a classical description, a system of kinks with a finite density 
forms a soliton lattice at zero temperature (see, e.g. [7]). However, the interaction between 
the solitons is short ranged and, at densities nS,o ^ 1, exponentially small. Quantum fluctu- 
ations, representing the Goldstone mode of the soliton lattice, destroy the periodic structure. 
Therefore, irrespective of the interaction constant, the transition to the Mott phase always 
occurs from a disordered phase, the latter described in terms of free spinless fermions. When 
the long-range Coulomb interaction is taken into account, the above described qualitative 
picture is modified. The unscreened repulsive Coulomb forces between the solitons lead to 
the formation of a Wigner crystal [8]. Strictly speaking, in a ID infinite chain of charges 
the 1/r-interaction does not fall off fast enough to remove the infrared divergence in the 
density-density correlation function. However, the Akp oscillating part of the correlator de- 
cays slower than any power [9], thus making it reasonable to speak of a nearly ordered state. 
At low densities nas ^ 1 {clb = ^ /e^rris and rus being the Bohr radius and soliton mass, 
respectively) the Wigner lattice with perod a = n~^ is not destroyed at exponentially large 
distances L < aexp{const/naB). Therefore, one can expect, in finite-size (mesoscopic) sam- 
ples, the "holons" to form a Wigner crystal, so that the transition to the insulating phase 
takes place from the Wigner crystal phase. We show that, in this case, the square-root 
singularity of the compressibility is changed to a logarithmic singularity. 

The holon density dependence of the behavior of the Wigner lattice is determined by 
the screening of the Coulomb interaction both at small and large distances. In quasi-one- 
dimensional systems, the ultraviolet cutoff is determined, as a rule, by the characteris- 
tic width A of a "one-mode" channel. Long-distance screening can occur when the one- 
dimensional chain of length L is situated close to a large metallic electrode, at at distance 
D <^ L. In this case the behavior of the compressibility k, = dn/dfi in the weak-screening 
{nD > 1) and strong-screening {nD <^ 1) regimes turns out to be qualitatively different. 
At weak screening, the dependence k = K,{n) does not differ from that in the case of an 
unscreened Coulomb interaction (with the obvious change L ^ D): K,{n) ~ l/ln(n-D). On 
lowering the density, one passes to the strong screening regime nD <^ 1, where the above 



logarithmic behavior is replaced by a power-law: K,{n) ~ {Dn)~'^. In the limit of arbitrarily 
small densities, this dependence would lead to the exponent 2/3 in the dependence of the 
compressibility on the chemical potential near the threshold: K^fi) r^ {fi ~ fic)~'^^^. How- 
ever, in the vicinity of the phase transition, quantum fluctuations melt the Wigner lattice 
and recover the universal square-root singularity K(fi) r^ [fi — fic)^^^'^, characterizing the 
commensurate-incommensurate transition [4]. The presence of a large metallic electrode 
allows the range of quantum fluctuations to be changed. On the other hand, it also leads to 
an additive renormalization of the critical value of //, /Xc = A — e^/AD, due to effects caused 
by image forces. 

The magnitude of the gap A in the spectrum of charged excitations in the insulating 
phase crucially depends on the character of the short-wavelength screening of the Coulomb 
interaction. This can be easily understood for weak Hubbard interaction U <^ t, when the 
charge excitations are quantum solitons with characteristic correlation length ^o ~ hcs/A (c^ 
is the velocity of the charge excitations in the Mott phase). If (^o/^)<^s ^ 1 (c^s = e^/ZiCg), 
the long-range Coulomb interaction weakly affects characteristics of the insulating phase, 
and the charge gap is determined by the well-known formula of Lieb and Wu [10]. 

When the opposite inequality is realized, {C,o/^)<^s ^ l, the Coulomb forces signiflcantly 
renormalize (increase) the Mott-Hubbard gap. It is physically evident that in this case the 
charge excitations — solitons — become heavier (A^ ^ A) due to the strong electrostatic 
energy, and "fragile" {^c ^ ^o) due to the unscreened Coulomb repulsion of charges "inside" 
the soliton. It seems reasonable to call these excitations Coulomb solitons, since their 
behavior in the conduction band at /i > A^; differs from that described above. Namely, 
at low densities, nC,.s <^ 1, the Coulomb solitons still condense into a Wigner crystal and 
display the logarithmic threshold singularity of the compressibility. However, at densities 
n^s ^ 1) but still much lower than a^^, the Wigner lattice is transformed to a sine-Gordon 
soliton lattice. The compressibility saturates at values k, ~ (Ac^s)^^, characteristic for a 
charge- density wave. 

Let us consider flrst the influence of long-range Coulomb forces on the properties of the 



metallic phase near the transition to the insulating state. In what follows, we shall assume 
that the transition is of the Mott-Hubbard type, although our conclusions remain valid for 
any one-dimensional charged system near the commensurate-incommensurate transition. 

Using the standard approach (see e.g. [11], and also [12]), one can treat the charge 
excitations of the half-filled, weak-coupling {U ^ t) Hubbard model as topological quantum 
solitons of the related sine-Gordon model (at /3^ = Svr), characterized by correlation length 
.^0 = hcs/A (here 2A is the Mott-Hubbard gap, and c^ — 2t + U/2tt is the velocity of 
charged excitations). The long-range repulsive forces will be taken into account by adding 
the screened (at short distances Ax ^ A) Coulomb interaction energy to the total energyof 
the system. The screening length A is determined, for example, by transverse dimensions of 
the confinement potential which forms the one-dimensional chain. 

Let A ^ .^0- Then it is reasonable to assume (the exact criterium will be given below) 
that the Coulomb interaction does not affect the intrinsic characteristics of the solitons {i.e. 
the Mott phase) but changes the character of interaction between them at large distances. 
Ax > A. Therefore, assuming the solitons to be point-like objects, we shall study their 
dynamics in the conduction band at low densities, n <^ a^^. At densities n < a^^ the 
Coulomb interaction leads to the formation of a soliton Wigner lattice with period a = n~^ 
{as = {hcs/e'^)C,o is the Bohr radius for charged particles with mass rris = A/c^). This 
ordered state is expected to be stable in chains of mesoscopically large length and therefore 
can significantly affect the properties of the metallic phase [13]. We shall study the influence 
of Wigner crystallization on the characteristics ofth e transition to the insulating (Mott) 
phase. 

Consider a situation, typical for "mesoscopic" experiments, when a massive metallic elec- 
trode is placed at the distance D {\ <^ D <^ L) from the chain. The electrode leads to 
screening of the Coulomb interaction at large distances. In this case, the classical (electro- 
static) energy of the lattice, with the electrode screening effect taken into account, equals 
[13] 
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where C is the Euler constant, and Ko{x) is the Macdonald function. The asymptotic 
behavior of the expression (1) in the hmits of high and low densities are 

Sw - eV ln{nD), nD > 1 (2) 

£w^n (-{^j + 2C(3)(eD) V, nD < 1. (3) 

Quantum corrections to the energy of the Wigner crystal (1) caused by zero-point fluctua- 
tions 

h r/"- ... 2 ^ (^'^^w 



are small, ASw <^ ^w (s is the sound velocity for the Wigner crystal). This is true at all 
densities n <^ a^^ in the weak screening regime, and within the interval a^D^^ <^ n <^ Oj^ 
at strong screening [nD <^ 1) as well. As could be expected, a strongly screened Coulomb 
interaction at low densities n < asD^"^ cannot stabilize the Wigner lattice which is already 
melted at zero temperature by zero-point fluctuations. In the absence of the macroscopic 
electrode, the regime Eq.(2) is realized, where D is replaced by the longitudinal length of 
the chain L. 

As already mentioned there is no true Wigner crystallization in an inflnite chain. 
Nonetheless, in mesoscopic-size systems. Coulomb correlations support almost perfect long- 
range order [9], making it possible to consider Wigner crystallization in the strict sense. 
Note that in the absence of impurities producing pinning of the Wigner crystal, the soliton 
lattice can freely slide along the one-dimensional channel. Therefore, in the ideal system, the 
response of a Wigner soliton-crystal to a low-frequency electric fleld is identical to that of a 
Fermi gas of solitons. (For instance, the charge stiffnesses are identical [14]). On the other 
hand, it is clear that the compressibilities of the two systems are different, so that, upon 
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crystallization of solitons, the nature of the singularity at the transition to the insulating 
phase on changing the chemical potential is changed. 

For the Fermi gas of solitons (with the rest energy A and mass rris = A/c^), the ther- 
modynamic potential of the system Q, calculated in the low-density limit, when the solitons 
can be treated as noninteracting particles, equals 

2 2 

f]p(n) = (A-/x)n + ^n3. (5) 

bnis 

Equation (5) leads immediately to the square-root singularity in the equilibrium density of 
holons near the phase transition to the insulating phase [5,6]: 

nifi) =0, /i < A (6) 



2 in I 

n{fi) = -^Jfi-A, fi>A. (7) 

TT Cs 

When Coulomb correlations are taken into account, the kinetic energy of the holons given 
by the last term in (5) should be changed by the electrostatic energy of the Wigner lattice 

(1). 

In the absence of the screening provided by the metallic electrode, the energy is given 
by Eq.(2), and instead of Eq.(5) we have for the thermodynamic potential 

f]^(n) = {A- fj,) n + e^n^ In(nL). (8) 

Comparing (8) with (5) and (7), we find that the square-root singularity of the compress- 
ibility is replaced by a logarithmic singularity: 

K-\fi)c^2e'\n ^^~^^^ . (9) 

We have dropped irrelevant numerical factors in the argument of the logarithm. It is assumed 
in Eq.(9) that the deviation Afi = /i — A is muchla rger than the Coulomb energy e^/L. One 
can easily check that, in this regime, the maximum length of the chain should not exceed 
Lm ~ aexp{const/naB) ^ a. For a larger length, L > L^, the ordered structure would be 
destroyed. 
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The presence of a massive metallic electrode makes the situation more complex. In 
the weak screening regime (2), similar to (9), we have a logarithmic dependence of the 
compressibility on the density 

K'\n) = 2e'^\n{enD), nD > 1. (10) 

On further decreasing the density, this dependence is replaced by a power-law, corresponding 
to the strong-screening regime: 

K-\n) = 24C (3) {enDf. (11) 

According to Eq.(ll), the dependence of the equilibrium density on chemical potential is 
given by 

n{fx) = 0, fi<A* (12) 

^- A* 



"(^) = k 





fl< /\^ 


1/3 






/i> A*, 



(13) 



where A* = A — e^ /AD. Substituting (13) into (11), one finds that the logarithmic singu- 
larity (9) of the compressibility, taking place for unscreened Coulomb potential, is changed 
to a power-law k^^ r^ (fi — A*)^/^. However, in the strong screening regime, at densi- 
ties n < aBD~'^, the Wigner lattice is destroyed locally by quantum fiuctuations. Therefore, 
Eq.(13) is valid at A/i = fi—A* > {e'^/D){aB/D)^. There is a region near threshold where the 
universal square-root singularity, characterizing the commensurate-incommensurate transi- 
tion, is recovered. Notice that the massive electrode leads to a shift of the critical value 
of the chemical potential, /ic = A — e^/AD, caused by the image forces (see e.g. [15]). 
This shift is small in the case of weak Coulomb coupling as = e'^/hcg ^ 1. In the strong- 
coupling regime, ag > 1, one should also account for strong Coulomb interactions affecting 
the insulating phase. 

For strong Coulomb coupling, as > I, the characteristic (Bohr) radius as turns out to be 
shorter than the correlation length ^o^ indicating the possibility of strong renormalization of 
the characteristics of single solitons by Coulomb interaction. From physical considerations 



one expects that, due to strong Coulomb interaction, solitons become heavier and extended. 
In other words, the effective correlation length C,s, being the size of the "Coulomb" soli- 
ton, may substantially esceed ^o- We shall study this situation assuming that screening is 
extremely small, A < .^o- 

A consistent solution of this problem suggests introducing the long-range Coulomb in- 
teraction at the level of the Hubbard model and calculation of the resulting spectrum of the 
charge excitations in such a system. However, for a half-filled system, repulsive Coulomb 
effects can only increase the charge gap. For this reason, reduction of the problem to the 
sine-Gordon model still remains reasonable. First we shall study the effect of Coulomb forces 
on topological solitons of the quasiclassical sine-Gordon model. Then we we shall present 
arguments explaining the applicability of the obtained results to the holon dynamics in the 
Hubbard model. 

In the presence of long-range Coulomb forces, the Lagrangian of the sine-Gordon model 
has the form: 



c = ^ 

Cs 



87r2 7-00 ^(x - yy + A2 ' 

We shall be interested in topological solitons (kinks) of the model (14). At j3'^ <^ 1 the 
model is quasiclassical, and the excitations we are interested in can be found by using trial 
functions. 

A trial function, describing a static topological soliton Aip = ip{x = +oo) — ip{x = 
— oo) = 27r//3 of the model Eq.(14), will be chosen in the form corresponding to the sine- 
Gordon model: 

Vcix) = - arctan f exp ( ^^ M j . (15) 

Here Xq is the center of the soliton, and d* is a variational parameter which determines the 
soliton size. Its value is found by minimizing the soliton rest energy: 
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Substituting (15) into (16) and doing the integrals, we easily find that, for a weak coupling 
as <^ 1, the parameter d* coincides with the size of an unperturbed kink, d* = do = Cs/wq, 
E{d*) = Eg = Shuo/P"^. For strong coupling, a^ > 1, the unscreened Coulomb interaction 
leads to multiplicative renormalizations of the size and energy of the soliton: 



d* = do [ a^ 



27r2 



E, 



c 



Shujo f l3^ 



P' 



'27r2 



lni/2 



do P' 



do_ 



P' 



— \las^^^ 



2tt' 



>rfo 



(17) 
(18) 



Let us show that expression (18) can be obtained within a consistent scheme, using solu- 
tions of the sine-Gordon equations. For a linearized problem, the Coulomb interaction can be 
taken into account exactly, leading to a modification of the spectrum of small perturbations 
("optical phonons") of the field ip (see, e.g. [9]) 



u\k) = uj', + k\l [ 1 + asf-^Ko{k\) ]. 



(19) 



Since at x > 1 -^'0(2^) is exponentially small, the Coulomb interaction mostly affects the 
long- wavelength dynamics, where the spectrum (19) takes the form 



u?{k) = cjg + Pc; 



«s/3' 



2„2 

' 27r2 



In 



k\' 



k\ < 1. 



(20) 



Expression (20) differs from the "phonon" spectrum of the unpertubed sine-Gordon 
equation in that the constant velocity c^ has been replaced by a momentum dependent 



effective "velocity" Cs(A;) ~ CsA/ln(l/|A;|A). Analysing the structure of the last term in 
Eq.(16), one concludes that, when studying the influence of weakly screened (A ^ rfg) 
Coulomb interaction on the topological soliton of the sine-Gordon model, it is sufficient to 
change the exact Lagrangian Eq.(14) by an approximate one in which the effects of Coulomb 
interaction are incorporated in a coordinate dependent velocity 
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~2^ "^ a" 



cl^c\x-Xo)=cl^ln^-:^—^. (21) 



In Eq.(21) it is assumed that \x — Xq\ ^ A. For this reason, it is possible to neglect in all 
calculations terms containing derivatives of the "velocity", \c'{x)/c{x)\ <^ 1. 

In the framework of such an "adiabatic" perturbation theory, the topological soliton has 
the standard form: 

ipc{x) = ±^ arctan[exp ( ^ ° J], (22) 

where the soliton size, do, is replaced by a coordinate dependent, smooth function d{x) = 
c{x)/ujq. One can easily check that the energy of the Coulomb soliton (22) exactly coincides 
with Eq.(18). 

Using the solution (22), the standard scheme of quasiclassical quantization can be easily 
developed (see, e.g. [16]). It can be shown that the one-loop quantum correction to the 
classical energy of the Coulomb soliton has the same form as in the usual sine-Gordon 
model, AEs = —hujo/ii. Within the traditional scheme of quasiclassical quantization of 
solitons, the relative smallness of quantum corrections is provided by the small value of 
the coupling constant /9^ ^ 1. In our case, solitons can be treated classically due to large 
Coulomb energy. Therefore, it seems natural to suggest that formula (18) remains valid in 
the extreme quantum case /3^ ^ 1 (from the point of view of the usual sine-Gordon model). 

When the charge sector of the half-filled, weakly repulsive {U/t <^ 1) Hubbard model is 
mapped onto the sine-Gordon model, the coupling constant /9^ turns out to be equal Svr [11]. 
The above analysis allows one to put forward a hypothesis that the long-range Coulomb 
interaction in the strong-coupling case e'^/hcs ^ 1 lead to a multiplicative renormalization 
of the Mott-Hubbard gap 



Notice that the appearance of fractional powers of a large logarithm is typical for various 
quantum problems that explicitly incorporate long-range Coulomb forces [13,9,17]. 
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Now we consider Coulomb solitons in the conduction band at /i > Ac. At low densities 
n^c ^ 1) where 

ec-eoay^ln^/^(|v^) (24) 

is the characteristic size of the Coulomb soliton, the charges in the conduction band can 
be considered as point-like. The unscreened Coulomb interaction leads to the formation of 
a Wigner crystal (see Eqs.(2),(9)). On increasing the density n > ^q^ <^ a^^, the solitons 
start to overlap strongly, forming a sine-Gordon lattice. 

The energy density of such a classical lattice can be readily estimated, using well-known 
periodic solutions of the sine-Gordon model 

,,,W = i{. + 2am||:^l}^ (25) 

Here am{z) is the elliptic amplitude, and k{n) is the elliptic modulus fixed by the soliton 
density 

2kK{k) = {nd*y\ (26) 

where K{k) is the complete elliptic integral of the /-order, and d* is defined in Eq.(17). The 
energy density of the soliton lattice equals (see also [18]) 

E,{n) = Ec^ {E{k) - i(l - e)K{k)} (27) 

{E{k) is the complete elliptic integral of the //-order, Ec being the energy of the Coulomb 
soliton). Expressions (25)- (27) are valid in our case at densities nd* > 1. In the limit of 
a dense soliton system nd* ^ 1 {k ^ 0), the soliton lattice smoothly transforms into a 
charge- density wave 

Epin) ~ ^iAsd*)n' (28) 

and the compressibility of the system is no longer dependent on the density: k^^ ~ A^ii* = 
const. 
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In conclusion we have shown that long-range Coulomb forces drastically modify the 
properties of ID electron systems in vicinity of the metal-insulator phase transition. In the 
metallic phase a weakly-screened Coulomb interaction leads to the formation of a Wigner 
crystal of charged quasiparticles and therefore affects the critical behaviour of the system 
at the transition point. The properties of the insulating phase are also changed if the short 
wavelength screening length is sufficiently small. In this case the Mott-Hubbard gap is 
strongly renormalized and the charged excitations in the Mott phase can be described as 
quasiclassical Coulomb solitons. 
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